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We investigate consequences of adding irrelevant (or less 
relevant) boundary operators to a (1 + l)-dimensional field 
theory, using the Ising and the boundary sine-Gordon model 
as examples. In the integrable case, irrelevant perturbations 
are shown to multiply reflection matrices by CDD factors: the 
low-energy behavior is not changed, while various high-energy 
behaviors are possible, including "roaming" RG trajectories. 
In the non-integrable case, a Monte Carlo study shows that 
the IR behavior is again generically unchanged, provided scal- 
ing variables are appropriately renormalized. 



In this paper we study the effect of adding irrelevant or 
less relevant scaling operators to a given theory, focusing 
on boundary perturbations in a (1 + l)-diinensional field 
theory. This is a problem of important physical interest. 
Consider for instance tunneling in the fractional quantum 
Hall effect, which has been intensely studied recently in 
the context of Luttinger liquids, shot noise and fractional 
charge measurements. For filling factor v — 1/3, there 
is a single relevant tunneling operator, but there are two 
of them for i' — 1/5. Most analytical calculations, in 
particular the exact solutions, do hold only when the 
second, less relevant operator has been scaled away. In 
practice, however, this operator will always be there, and 
it is important to be able to evaluate its role. 

The standard expectation is simply that one can ne- 
glect all less relevant and irrelevant operators when, com- 
puting the properties close to the IR fixed point .El This 
is, however, entirely based on w!ak coupling expansions, 
which can be quite misleading.! For instance, the added 
perturbations could become relevant at strong coupling, 
entirely changing the physics in a non-perturbative way. 

Our study supports the standard expectation in a va- 
riety of cases related with the tunneling problem. On the 
conceptual side, one of our main results is that, in the 
integrable case, additional irrelevant perturbations con- 
tribute so-called CDD (Castillejo-Dalitz-Dyson) factors 
to the boundary reflection matrices R. This is not so 
surprising: as discussed in Ref. ||, R is fully determined 
up to such CDD factors by the boundary analogues of the 
Yang-Baxter equation, unitarity and crossing symmetry, 
and our results thus provide a natural explanation for the 
CDD ambiguity. Nevertheless, unexpected features are 
found close to the UV fixed point which can be altered 
by the added perturbations and is approached along a 
"roaming" renormalization group (RG) trajectory that 
oscillates between different boundary conditions. The 



relation between irrelevant perturbations and CDD fac- 
tors has been discussed independently for the bulk case 
in Ref. § 

To start, we consider the simplest possible situation, 
which is the scaling limit of the Ising model with zero ex- 
ternal field but with a boundary magnetic field h. This 
integrable napdel is described by a free Majorana fermion 
field theoryE Since the boundary spin operator coincides 
with the free fermion operator (in sharp contrast to the 
bulk case), the problem can be handled by solving the 
boundary equation of motion, which in turn determines 
the boundary reflection matrices. It is then easy to in- 
clude additional irrelevant operators: all perturbations 
that are quadratic in the fcrmions will still lead to a solv- 
able model. 

We illustrate this for a boundary perturbation made 
up of two terms, the spin operator and the stress energy 
tensor T = n : ipdip :, where d = dz and z = x + iy. The 
Euclidean action is 



A = 



dy 



dx (^tpdip — tpdip -\- miptp^ 
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where a is a fermionic boundary degree of freedom. The 
case A = was studied in Ref |[ When h = X = 0, the 
model is at the fixed point corresponding to free bound- 
ary conditions, or in the fermion language, — ip. It 
follows that on the boundary, T — —in : ipdy^ :, and the 
last term in Eq. (0) reads (A/tt) Jdy{T-\-T). 

Let us now determine R. The equations of motion at 
the boundary are 



(2) 



The bulk theory in Eq. (|i|) contains one type of particle, 
namely the free fermion A with mass m. The fermion 
fields can be expressed in terms of the particle-creation 
operator A'l'(0),H 
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where w = exp(i7r/4), = exp(— i7r/4). Inserting this 
expansion into Eq. and using the defining relation 
A\6) = R{9)A\-9), it is straightforward to obtain the 
reflection matrix. For simplicity, we focus on the mass- 
less limit, m — > 0. In that limit, we boost the rapidities 9 
such that the energy of the particles m cosh 61 exp(/?), 
while the momentum m|Sinh0 — > ±exp(/?) for right- or 
left-movers, respectively.u Physical rapidities are charac- 
terized by exp(/3) « T, where T denotes the temperature. 
The reflection matrix for rapidity (3 then reads 



R([3) 



This can be recast into the form 
R(P) = ~i tanh 



^~/^-^]cothf^~^^ 
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where (the -|- sign corresponds to f3g) 



exp(/3^'^) 



Vl + 4A/i2 ± 1 
4A 
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Sending the additional irrelevant coupling to zero, A — > 
0, we obtain exp(/3^) « while exp(/3|) w 1/2A 

diverges. Hence only the first term of the R matrix (Q) 
matters, while the second factor saturates at minus one: 
the boundary magnetic field h induces a fiow from free to 
fixed boundary conditionsp or, in the fcrmion language, 
from ip = ip to ip — Close to the UV fixed point 

(/3 — > oo), the reflection matrix is Rfree = i- Close to the 
IR fixed point (/3 — oo), we instead get Rfixed = —i- 

Let us now discuss the meaning of Eq. (^) . In the limit 
where both bare couplings are small, we have exp(/3|j) <JC 
exp(/3^). The vicinity of the free fixed point now corre- 
sponds to temperatures exp(/3]j) ^ T ^ exp(/3^). In 
that regime, the coth-factor in Eq. is close to mi- 
nus one, and Eq. (|l|) describes an Ising model with free 
boundary conditions perturbed by a relevant (cx h) and 
an irrelevant (cx A) term. There are now two possibilities. 
If we look at higher energies, i.e., run the RG backwards, 
the contribution of the irrelevant term is seen to grow 
while the the relevant term disappears. Equation (H) for 
the R matrix shows that in this limit the model flows to 
fixed boundary conditions, R — > —i. Therefore the UV 
fixed point is not characterized by free boundary condi- 
tions anymore, as one might have expected. If instead 
we look at lower energies, the irrelevant term is scaled 
away while the relevant term grows. From Eq. (|^) we see 
that the model again fiows to fixed boundary conditions, 
R —i. The result of running the RG forward (to lower 
energies) is very natural and confirms the general idea 



that irrelevant couplings should scale away. The result 
of running the RG backwards (to higher energies) is less 
natural, and its simplicity is probably connected to the 
fact that we are dealing with a free theory. In more com- 
plex cases, we do not expect that such a backward RG 
trajectory ever converges. 

The structure of Eq. (^) generalizes to more complex 
mixtures of N perturbations that preserve integrability. 
Up to a factor —i, the R matrix should thereby result as 
a product of N (minus) hyperbolic tangents with asso- 
ciated characteristic "Kondo" temperatures, exp{/3g) ^ 
exp(/?^) ^ . . . <C exp(/3^). If one considers some phys- 
ical property as a function of temperature, it will ex- 
hibit "non-universal" properties for T >, exp(/3^). For 
T <^ exp(/9g), however, all the hyperbolic tangents but 
the first one are saturated, and then the refiection matrix 

becomes R ~ ztanh {^^—^ ~ T")' identical to the case 
of a single relevant perturbation. Note that in these cal- 
culations, no cut-off has ever appeared. This is because 
of the implicit dimensional (normal-order) regularization 
used in integrable systems. In practice, the theory must 
often be regularized using a UV cut-off, whence a differ- 
ent renormalization of by the irrelevant couplings will 
occur. 

Running the RG backwards from physical tempera- 
tures in the range exp(/3]j) <C T <C exp(/3^) gives in- 
triguing results. If A'' is even, the model starts at free 
boundary conditions and ends up at fixed ones. The UV 
fixed point is approached by "roaming", where the RG 
trajectory comes again N/2—1 times close to free bound- 
ary conditions. For instance, if = 4, the backward RG 
trajectory reads: fixed <— free ^ fixed ^ free. If N 
is odd, the model ends up at free boundary conditions, 
and comes close to fixed boundary conditions {N — l)/2 
times. For instance, if iV = 3, the backward trajectory 
reads: free <— fixed <— free. Of course, in a given realistic 
situation, we expect the equivalent of N to be infinite, 
so the backwards trajectory will keep roaming forever as 
energies get higher and higher. In a cut-off regularized 
theory, such a behavior will only be meaningful up to 
energy scales given by the UV cut-off. In contrast, the 
forward trajectory always describes the same flow from 
free to fixed boundary conditions. 

The above considerations for the Ising model can now 
be applied to the boundary sine-Gordon (BSC) model, 
which in turn is related to the problem of a quantum 
impurity in a Luttinger liquid,l3 

A^- dx dy [(a,<i>)2 + (a,$)2] (6) 
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Here $(x,?/) is a folded bosonic field.0 Remarkably, the 
g — 1/2 BSG model can be mapped on a pair of 
Ising models.l3'B The first of these Ising models has fixed 
boundary conditions [h = oo), while the other sees an ef- 
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fective boundary magnetic field proportional to the cou- 
pling h. This mappins-|Can be made explicit at the level 
of reflection matricesB At g — 1/2, the BSG model 
has kinks and antikinks with a bulk scattering matrix 
S = —1. At the boundary, they elastically scatter one 
by one without particle production. The reflection am- 
plitudes for kink- kink scattering, J-" = R++, and kink- 
antikink scattering, Q = , areO 

P + Q = ~i, P - Q = z tanh (^^^^ - , (7) 

where exp(/3B) cx h^. The mapping onto the Ising mod- 
els is obtained by taking symmetric and antisymmetric 
linear combinations of kink and antikink, which then cor- 
respond to the fermions tp and ^ of the two Ising models. 
Clearly, the reflection matrices are consistent with 
the Ising R matrix (jj). 

Suppose now that we add a general combination of 
the two irrelevant operators cos[2-y/27r(7$(x — 0)] and 
: [dy'^ix = 0)f :. At 5 = 1/2 and h ^ 0, both op- 
erators are irrelevant with scaling dimension 2. A gen- 
eral combination of this sort can now be rewritten in 
terms of the underlying two Ising models. By means 
of bosonizatiora we obtain the equivalent combination 
A : ipdyip : +fj, : £,dy^ :. By using the same reasoning as 
in the previous paragraph, Eq. (0) yields 



P + Q = i coth 



P — Q = —i tanh 
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In the limit of small bare couplings, exp(/3]j) cx 
/i^,exp(/3^) oc l/A,exp(/3^) oc l/fi, and therefore ^ 
(3%,P%. For temperatures T < exp(/3^), exp(/3|), the 
coth-f actors are saturated, and we recover Eq. (0). For 
high temperatures, in contrast, one gets P + Q i, 
P — Q ^ —i, corresponding to Neumann (free) bound- 
ary conditions. The situation (||) is thus similar to the 
case = 3 for the Ising model. Notice however that in 
the special case ^ = 0, we always get P + Q = —i, and 
the situation now looks like the case A^ = 2 for the Ising 
model, with identical Dirichlet (fixed) boundary condi- 
tions in the UV and IR. The Neumann boundary con- 
dition is generally found only in the temperature regime 
exp(/?]3)«r«exp(/3|). 

Following Ref. the exact transport properties can 
be computed for arbitrary boundary interactions once 
the reflection matrices are knows. The 5 = 1/2 linear 
conductance in units of ge^/h isllll 



df3- 



IQ(/3 + inr)|2 



(9) 



where Q is given by Eq. (||). This formula results from 
Ref. |o| after an integration by parts, and holds provided 
IQP is not a constant, so there is a flow indeed. Figure |l| 



shows plots for Eq. (||). In the absence of the irrelevant 
couplings, the conductance goes from the perfect value 
ge^ /h at high T (Neumann boundary conditions) down 
to G = as T ^ (Dirichlet). With the irrelevant 
couplings, the situation is more complex. For /i = 0, in 
the UV limit Dirichlet boundary conditions [G = 0) are 
approached again. However, for finite /i, the UV fixed 
point corresponds to Neumann boundary conditions, see 

Fig. g. 
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FIG. 1. Temperature dependence of the (7 = 1/2 conduc- 
tance for h = 1 and various A, /j,. 

Let us now discuss the effect of irrelevant operators in 
a cut-off regularized theory. We first consider the g — 1/2 
BSG model (^) and only add the first of the above two 
operators, cos[2y/2TTg^{x — 0)], with a coupling constant 
A. Its effect is studied by Monte Carlo simulation pfor the 
pinning function F{x) of the Friedel oscillation.E3 The 
2A:ir-oscillatory charge screening cloud p{x) of a spinless 
Luttinger liquid around the impurity at position a; = is 



Pix) - po 



cos{2kpx + riF)(a/xy F{x) , 



(10) 



where po = kp/n and 77^ is a phase shift. In the ab- 
sence of the irrelevant coupling, F depends on x/xb 
only, where xb = a{ha)~'^/^^~^^ is a lengthscale corre- 
sponding to exp(/3B). For T = 0, energy scales are re- 
lated to the distance x with x ^ xb meaning low ener- 
gies. In Eq. (|l^) we have introduced the lattice spacing a 
which, sai'ves as a high-energy cut-off. General arguments 
showB'EHl that F[x) ^ 1 as x ^ 00. Figure || shows nu- 
merical results for the g =\/2 pinning function£3 From 
the A ^ data, one sees that in the low-energy portion 
X ^ Xb , the additional irrelevant operator can simply be 
incorporated by a rescahng of the lengthscale Xb ov the 
associated rapidity (Sb (such a rescaling could also have 
been used in the fits of figure 1). A good estimate for 
the observed size of the rescaling can be pbtained from 
a self-consistent harmonic approximation.E3 At high en- 
ergy scales, X <^xb^ however, the curves do not scale on 
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top of each other, in agreement with the general picture 
outhned above. 
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FIG. 2. Monte Carlo results for the g = 1/2 pinning func- 
tion {h = 1). A hard energy cut-off ojc = 1/a was employed, 
and the simulations were done at very loW|but finite tempera- 
ture LJc/T = 1200. The exact T = resultO for A = is given 
as dotted curve. Data for A = are given as solid curve, and 
for A = 0.6 as dashed curve after rescaling a;s by a factor 
0.43. Statistical errors are of the order of the thickness of the 
curves. 



to resort to numerical simulations. We have studied the 
combination of the usual BSG term cos \^/2¥g<^{x = 0)] 
and of cos[2-/27r(7<i>(2; = 0)], which is not integrable away 
from g = 1/2. In Fig. ^, we show data for g = 1/5, where 
the second operator is relevant (with scaling dimension 
4/5), albeit of course less relevant than the first. We ob- 
serve that at low energy scales (x ^ xb), the coupling 
A has again no effect besides a renormalization of xb- 
The scaling form of F is altered only for x ^ xb, giving 
rise to a picture that appears to be quite similar to the 
integrable one. 

Our results should not lead to the belief that irrele- 
vant operators are always of little importance. The in- 
tegrable situation is quite special of course, since it does 
not even require the introduction of a cut-off. In non- 
integrable cases, the necessary cut-off might give rise to 
more complex behaviors. More importantly, the "roam- 
ing" RG trajectories we observed when approaching the 
UV fixed point in the Ising model should serve as a warn- 
ing that there are situations where irrelevant operators 
lead to counter-intuitive behavior, sometimes even chang- 
ing the nature of the IR fixed point. In the boundary 
sine-Gordon model, however, this seems never to be the 
case. 
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FIG. 3. Monte Carlo results for the g = 1/5 pinning 
function (h — 1). Data for A = are given as the solid curve, 
and for A — 0.6 as the dashed curve after rescaling x_b by a 
factor 0.19. 

The inclusion of additional irrelevant or less relevant 
terms is not restricted to g = 1/2. First, recall that a 
general consequence of integrability is the existence of a 
whole infinity of quantities in involution. A model per- 
turbed by a combination of all these quantities is then 
also integrable, and the associated R matrices should dif- 
fer from the standard ones only by the inclusion of CDD 
factors, giving rise to results similar to the Ising ones. As 
for adding operators that destroy integrability, we have 
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